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Abstract 

I discuss some recent developments in chiral effective field theory for few-nucleon systems. 
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1. Introduction 

Almost two decades ago Weinberg formulated his ideas to extend chiral perturbation 
theory to systems with several nucleons [1]. This seminal work initiated extensive re- 
search towards better understanding of the structure of nuclear forces and the dynamics 
of few- and many-nucleon systems based on the powerful framework of effective field the- 
ory (EFT). In this talk I outline the current status of chiral EFT for nuclear forces and 
applications in the few-nucleon sector. I discuss the structure of the subleading contribu- 
tions to the three-nucleon force (3NF) and the leading four-nucleon force (4NF) which 
are currently being investigated. The role of A-excitation will also be addressed. The last 
part of this talk concerns lattice simulations of light nuclei. 

2. Chiral EFT for few nucleons: present status 



Chiral EFT is an appropriate framework to analyze the properties of few-nucleon 
systems at low energy! 1 1 It is based on the most general effective Lagrangian for Goldstone 
bosons (pions) and matter fields (nucleons and perhaps deltas) consistent with the chiral 
symmetry of QCD. For energies below the pion-production threshold, a nonrelativistic 

Email address: e. epelbaumafz-juelich.de (E. Epelbaum). 
1 At very low energies, it is advantageous to use pionless EFT, see [2] for a recent review article. 
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treatment of the nucleons is justified. Further, it is advantageous to eliminate the pionic 
degrees of freedom which gives rise to the (nonlocal) nuclear Hamiltonian 

ffnucl = H + V 2N + V 3N + V 4N + ... . (1) 

The resulting quantum-mechanical A-body problem can be dealt with using rigorous 
few-body techniques or many-body methods. 

The derivation of the nuclear potentials from the effective Lagrangian can be carried 
out straightforwardly utilizing the so-called chiral expansion. The importance of a partic- 
ular contribution to the nuclear Hamiltonian is determined by the corresponding power 
of the expansion parameter Q/A where Q and A refer to the generic low- momentum 
scale associated with external nucleon momenta or Mj and the pertinent hard scale, 
respectively [1]. For the two- nucleon force (2NF), this expansion has the form 

V2N - V 2N + V 2N + V 2N + V 2N + . . . , (2) 

where the subscripts refer to the chiral order. The leading contribution V 2 ^ is due to Itt- 
exchange (OPE) and two NN contact interactions without derivatives. The corrections 
result from contact interactions with increasing number of derivatives and/or insertions 
of M% and exchange of pions. In particular, 27r-exchange starts to contribute at next- 
to-leading order (NLO) and provides an important ingredient of the nuclear force [3]. 
At present, the 2NF has been worked out and applied in the NN system at next-to- 
next-to-next-to-leading order (N 3 LO) [4,5], see [6] for the pioneering work along this 
line. Its long-range part involves contributions from In- 2-k- and 37r-exchange and is 
parameter-free since the corresponding low-energy constants (LECs) are known from nN 
scattering. Numerically, the 37r-exchange potential turns out to be negligibly small [7]. 
The short-range part consists of 24 independent contact interaction^ whose strengths 
were determined from fits to NN low-energy data, see Ref. [8] for the discussion on 
resonance saturation of the corresponding LECs. Both N 3 LO calculations of Ref. [4] and 
[5] yield compatible results for the NN system (within the theoretical uncertainty) and 
demonstrate an accurate description of the low-energy scattering data and the deuteron 
properties. 

Systems with three- and more nucleons provide an excellent testing ground for chiral 
nuclear forces and were extensively studied in this framework over the past few years. 
Chiral power counting explains naturally the observed hierarchy of nuclear forces V2N ^ 
V 3 n > ViN ■ ■ ■ with 

V 3N = V$+V$ + ..., V iN = V^ + V^ + ..., .... (3) 

The first nonvanishing contribution to the 3NF appears at N 2 LO relative to the leading 
2NF and results from 27r-exchange, l7r-exchange with 2N contact interaction and the 
purely short-range 3N contact term. While the 27r-exchange contribution is parameter- 
free, the two other topologies depend on the unknown LECs D and E which were fitted 
to the 3 H binding energy and nd doublet scattering length [9] and to the combinations of 
the 3 H/ 4 He [10] and 3 H/ 10 B [11] binding energies. The resulting nuclear Hamiltonian was 
extensively tested in the 3N continuum [12,13,14,15] and applied to study the properties 

2 This number does not include isospin-breaking terms and refers to Weinberg's counting rules based 
on the naive dimensional analysis. Alternative counting schemes are currently being explored. 
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Fig. 1. Differential cross section (in mb/sr) for elastic Nd scattering at 10 MeV (left panel) and 65 MeV 
(middle panel) at NLO (light-shaded bands) and N 2 LO (dark-shaded bands). References to data can be 
found in [9]. Right panel shows pd breakup cross section (in mb MeV -1 sr~ 2 ) along the kinematical locus 
S at 19 MeV nucleon energy in the SCRE configuration with a = 56° [12]. Dashed and dashed-dotted 
lines are results based on the CD Bonn 2000 2NF [17] combined with the TM99 3NF [18] and the coupled 
channel calculation including the explicit A and the Coulomb interaction [19], respectively. 

of light nuclei using the no-core shell model approach [10,11,16]. The calculated 3N scat- 
tering observables are generally in a reasonable agreement with the data which improves 
when going from NLO to N 2 LO, see e.g. the left and middle panels in Fig. 1. Notice 
however that the theoretical uncertainty at N 2 LO appears to be rather large already at 
moderate energies, see the middle panel in Fig. 1. We further emphasize that there are 
certain observables for which remarkably large discrepancies with the data are found even 
at very low energy. A well-known example is the so-called symmetric space-star (SST) 
Nd breakup configuration at En = 13 MeV [20,9], in which the plane in the cm. system 
spanned by the outgoing nucleons is perpendicular to the beam axis. Recently, pd data 
for a similar symmetric constant relative-energy configuration have been measured at 
Ed = 19 MeV [12]. This geometry is characterized by the angle a between the beam 
axis and the plane in the cm. system spanned by the outgoing nucleons. Similar to the 
SST geometry, one observes large deviations between the theory and the data, see right 
panel in Fig. 10 which also hold for calculations based on modern phenomenological po- 
tentials. The included 3NFs have little effect on the cross section while the effect of the 
Coulomb interaction is significant but removes only a part of the discrepancy. For more 
details on these and related topics the reader is referred to Ref. [21], see also Rcf. [22]. 

3. N 3 LO contributions to three- and four-nucleon forces 

Although the results for few-nucleon observables already look quite promising at N 2 LO, 
it is mandatory to extend these studies to N 3 LO as it has already been done for the 2N 
system. This would allow one to test the convergence of the chiral expansion and might 
shed some light on the long-standing puzzles in the 3N continuum. The extension to 



3 We emphasize that chiral EFT results are obtained without taking into account Coulomb interaction. 
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Fig. 2. Different 3NF contributions at N 3 LO. Solid and dashed lines represent nucleons and pions, 
respectively. Solid dots denote the leading-order vertices. 

N 3 LO requires the incorporation of the leading corrections to the 3NF which feed into 
five different topologies, see Fig. 2, and are currently being worked out. Interestingly, 
there are no contributions from higher-order vertices in the effective Lagrangian (except 
I/m-corrections) so that the 3NF at N 3 LO is parameter- free. Notice further that the 
spin-space structure of the chiral 3NF at N 3 LO is much richer than the one of the 
currently available 3NF models. Partial results are already available for certain types of 
contributions. In particular, one-loop and relativistic 1/m-corrections to the 27r-exchange 
topology were considered in Refs. [23] and [24], respectively. The results of Ref. [23] based 
on the AV18 2NF [25] indicate that the one-loop 27r-exchange corrections have little effect 
on 3N scattering observables. 

In addition to the 3NF corrections, four-nucleon force (4NF) starts to contribute at 
N 3 LO in the chiral expansion. It has been worked out recently using the method of 
unitary transformation [26]. The leading 4NF results from diagrams depicted in Fig. 3 
where only graphs are shown which yield nonvanishing contributions. It is governed by 
the exchange of pions and the lowest-order nucleon-nucleon contact interaction and in- 
cludes effects due to the nonlinear pion-nuclcon couplings and the pion self-interactions 
constrained by the chiral symmetry of QCD. The obtained 4NF is local and does not 
contain any unknown parameters. The expectation values of the individual 4NF contri- 
butions obtained in the pioneering study of Ref. [27] are typically of the order of few 
hundreds of keV, which agrees with estimations based on dimensional arguments. One 
also observes strong cancellations between individual contributions. 
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Fig. 3. 4N diagrams generating nonvanishing 4NF contributions at N 3 LO. Diagrams which result from 
the interchange of the nucleon lines and/or application of the time reversal operation are not shown. For 
remaining notation see Fig. 2. 

4. Role of A-excitations 




The A-isobar is known to play an important role in hadronic and nuclear physics due 
to its low excitation energy Am = wa — m = 293 MeV, and strong coupling to the 
ttN system. It can, therefore, be argued that the explicit inclusion of the delta allows 
one to resum a certain class of important contributions and thus leads to an improved 
convergence as compared to the delta-less theory. Indeed, it is well known that e.g. the 
LECs c 3i 4 accompanying the subleading irirNN vertices are saturated by the A-isobar 
which leads to their rather large numerical values [28] . As a consequence, the subleading 
27r-exchange 2NF turns out to be much stronger than the leading one [29]. A similar 
tendency was also found for the 37r-exchange [7] and the charge-symmetry breaking 2ir- 
exchange potentials [30]. In EFT with explicit deltas, the dominant portion of these large 
subleading contributions is shifted to the lower chiral order leading to a more natural 
convergence pattern. 

The explicit inclusion of the delta in chiral EFT can be achieved via the so-called small 
scale expansion (SSE) [31] in which Am is treated as an additional small parameter. The 
leading contributions to the 2NF due to intermediate A excitations arise at NLO from 
diagrams shown in the first raw of Fig. 4 and were considered in [6,32]. In [33] we worked 
out the subleading contributions which are generated at N 2 LO by diagrams depicted 
in the second raw of Fig. 4. To obtain numerical results, one needs to determine the 
combination of the the subleading irNA LECs 63 + b% which enters at this order and 
the LECs Cj whoes values differ from the ones obtained in the delta-less theory. This 
was achieved by fitting the threshold coefficients of ttN scattering calculated at second 



Fig. 4. NLO (first raw) and N 2 LO (second raw) contributions to the 2-7r-exchangc 2NF with single and 
double A excitations. Double lines represent A-isobars, filled squares denote subleading vertices. For 
remaining notation see Fig. 2. 
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Fig. 5. Isoscalar central (Vc(r)) and isovector spin-spin (W,$(r)) potentials in r-space in EFT with 
(left panel) and without (middle panel) explicit A's. Right panel: 3 i<3 NN phase shift in EFT with and 
without explicit A's calculated in first Born approximation. The filled circles (open triangles) depict the 
results from the Nijmegen multi-energy PWA [34] (Virginia Tech single-energy PWA [35]). In all cases, 
spectral-function regularization with the cutoff A = 700 MeV has been used. 



order in the SSE. As expected, we found that the values of the LECs C34 are reduced 
in magnitude in the theory with explicit A's. Using the SU(4) (or large N c ) value for 
the leading irNA LEC Ha = 3<7a/(2-\/2) = 1.34 with gA being the nucleon axial-vector 
coupling, we found C3 = —0.79 GeV^ 1 and C4 = 1.33 GeV -1 which has to be compared 
with c 3 = -3.87 GeV" 1 and c 4 = 2.89 GeV" 1 in the delta-less EFT. 

In Fig. 5 we compare the isoscalar central and isovector spin-spin 27r-exchange contri- 
butions at NLO and N 2 LO obtained in EFT with and without explicit A's. Contrary to 
the delta-less theory where the whole contribution is generated at N 2 LO by subleading 
27r-exchange, the dominant contribution in the delta-full theory results at NLO from the 
leading 27r-exchange graphs. This much natural convergence pattern is also observed for 
peripheral phase shifts, see Fig. 5 for an example. 



5. Few nucleons on a lattice 



The last topic I would like to address concerns nuclear lattice simulations using chiral 
EFT [36], see also [37] for a similar work based on pionlcss EFT. In this framework, the 
effective Lagrangian is formulated on a Euclidean lattice and the path integral is evalu- 
ated by Monte Carlo sampling. Pions and nucleons are treated as point-like particles on 
the lattice sites, and n times the inverse lattice spacing sets the cutoff scale in momentum 
space. By using hadronic degrees of freedom and concentrating on low-energy physics, 
it is possible to probe larger volumes, lower temperatures, and far larger numbers of 
nucleons than in lattice QCD. In [36] simulations were carried out at lowest order in 
the chiral expansion. The sign problem is suppressed due to the SU(4) positivity of the 
lattice path integral for any even number of nucleons and the approximate SU(4) sym- 
metry of the low-energy nuclear interactions. Having fixed three unknown parameters 
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Fig. 6. Left panel: antisymmetric linear combination of the nucleon density correlations for the deuteron 
in the yz-plane. The deuteron spin points in the +z-direction. Right panel: CPU time versus the number 
of nucleons, A, measured relative to the A = 2 deuteron system. 

entering the 2NF using the 2 H binding energy, the NN 1 Sq scattering length and the 
average S-wave effective rangeQ we computed various properties of light nuclei, see left 
panel in Fig. 6 for one example. The results for the deuteron properties and the obtained 
value for the triton binding energy Ea^ = —8.9(2) MeV are within 5% of the data while 
the a-particle binding energy Ei He = —21.5(2) MeV is about 25% smaller in magnitude 
than the experimental value of —28.296 MeV. Encouraging results were also obtained for 
the computational scaling showing for A < 8 an approximately linear dependence on A, 
see right panel in Fig. 6. 

6. Summary and outlook 

In this talk I discussed the structure of the nuclear force in few lowest orders in the 
chiral expansion. In the 2N system, accurate description of the low-energy data is achieved 
at N 3 LO. The results for 3N scattering and the properties of light nuclei are so far only 
available at N 2 LO. While most of the calculated 3N scattering observables are in a 
reasonable agreement with the data, the theoretical uncertainty at this order appears to 
be rather large. It is therefore mandatory to extend this studies to N 3 LO. The corrections 
to the 3NF which arise at this order are currently being worked out. The leading 4NF 
which also contributes at N 3 LO has already been derived. 

In the second part of my talk I discussed the 2NF due to intermediate A-excitations at 
N 2 LO in the SSE. It is demonstrated that the explicit inclusion of the A-isobar in chiral 
EFT allows to improve the convergence of the low-momentum expansion. In the future, 
this work should be extended to study the role of the A-isobar in the 3NF and isospin- 

4 The additional third parameter corresponds to higher-order contact interactions which were included 
to overcome a clustering instability, see [36] for more details. 
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breaking interactions. It would also be interesting to work out N 3 LO contributions in the 
SSE. 

Finally, in the last part of my talk, I discussed the leading-order results for the prop- 
erties of light nuclei based on the lattice formulation of chiral EFT. A generalization 
to higher orders in the chiral expansion and to scattering observables, see [38] for the 
important step in this direction, is in progress. 
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